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The binary companion mass ratio distribution: an imprint 
of the star formation process? 
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ABSTRACT 

We explore the effects of dynamical evolution in dense clusters on the companion 
mass ratio distribution (CMRD) of binary stars. Binary systems are destroyed by 
interactions with other stars in the cluster, lowering the total binary fraction and 
significantly altering the initial semi-major axis distribution. However, the shape of the 
CMRD is unaffected by dynamics; an equal number of systems with high mass ratios 
are destroyed compared to systems with low mass ratios. We might expect a weak 
dependence of the survivability of a binary on its mass ratio because its binding energy 
is proportional to both the primary and secondary mass components of the system. 
However, binaries are broken up by interactions in which the perturbing star has a 
significantly higher energy (by a factor of > 10, depending on the particular binary 
properties) than the binding energy of the binary, or through multiple interactions in 
the cluster. We therefore suggest that the shape of the observed binary CMRD is an 
outcome of the star formation process, and should be measured in preference to the 
distributions of orbital parameters, such as the semi-major axis distribution. 

Key words: stars: formation - kinematics and dynamics - binaries: general - open 
clusters and associations: general - methods: numerical 



1 INTRODUCTION 

' One of the major unsolved questions in astrophysics is 
whether the star formation process is universal, or whether 
it varies as a function of the local environment. Stud- 
ies of the Initial Mass Function (IMF) in nearby star 
forming regions have shown that the IMF is invariant 
l|Bastian, Covey fc Mever1l2010l ). 

A further exploration of this universality hypothesis 
can be undertaken by studying t he properties of binary 
stars in different e nvironments (e.g. iReggiani fc Meverll201ll ; 
King et al.ll2012a .b"). If the properties of binaries vary as a 
function of environment, then we can conclude that star 
formation is not universal, and depends on some localised 
parameter. 

However, comparing binary statistics in different star 
formation environments is far from straightforward. Firstly, 
one needs to compare the statistics of binaries over a com- 
mon primary mass range, separat ion range and mass ratio 
sensitivity (e.g. iKing et al.ll2012al ). 

Secondly, dynamical evolution in dense star formation 
regions has been shown to alter the overall binary fraction, 
and the separation a nd orbital eccentricity distributions 
(e.g. lKroupalfl995al rbl: iKroupa. Petr fc McCaughreanlll999l ; 
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Parker. Goodwin fc Allison! l201ll : iMarks. Kroupa fc Ohl 



20111 ). It is possible to 'reverse engineer' the initial condi- 



tions of star formation in clusters to determine the amount 
of dynamical processing that has taken place, but in order 
to do this one must first make assumptions about the 
initial bin ary properties, and the highes t dens it y the cluster 
attai ned (|Kroupal Il995al : iParker et all 120091 : iMarks et all 
I2011T ). Such assumptions are non-trivial if the dynamical 
processing has occurred ' locally' (e.g. withi n pockets of 
substructure in the cluster. IParker et al.ll201ll ). rather than 
globally. 

A more fruitful approach would be to focus observa- 
tional effort on a measurable parameter of binary stars that 
may not be affected by dynamical evolution. One such pa- 
rameter could be the distribution of binary mass ratios, 
q, where q = m s /m p , and m p and m B are the masses of 
the primary and secondary components of th e system, re - 
spectively. Recent work by IReggiani fc Meyer! l|201ll . 120131 ) 
compared this companion mass ratio distribution (CMRD) 
in different star formation regions, and open clusters, with 
binaries in the Galactic field. They find that the shape of 
the CMRD is consistent with being linearly flat in the field 
(dN/dq oc g-°' 50 ±°' 29 ) ] anc j m mos t star forming regions 
(with the possible exception of the Taurus association). 

Dynamical interactions will reduce the overall fraction 
of binary systems in a cluster, but it is unclear whether 
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the shape of the CMRD is altered. In iV-body simula- 
tions of_^ery_l_ow_jriass binari es (VLMBs) in dense clus- 
ters, IParker fc Goodwinl (|201ll ) found that the CMRD for 
VLMBs is uniformly lowered, with no preference for destroy- 
ing systems with a low q. However, VLMBs typically form 
as equal-mass systems, and it is unclear how the CMRD 
would change for binaries with higher primary masses and 
an initially flat CMRD, as is observed in the field. 

In this paper, we investigate the effects of dynamical 
interactions in dense clusters on the shape of the CMRD. We 
summarise the analytics of binary destruction in Section [2] 
we outline our numerical simulations in Section|3] we present 
our results in Section 3] we provide a discussion in Section [5] 
and we conclude in Section [6] 



The outcome of an interaction between a single star 
and a binary depends on the relative velocity of the per- 
turbing star, upcrt, th e critical velocity, v c , and the or- 
bital velocity, n or b (e.g. iFregeau et al.ll2006l . and references 
therein). If vp crt < v c , then destruction is not possible, 
but such a scenario favours exchange interactions (which 
would likely increase the mass ratio). On the other hand, if 
V c < "Pert < forb, then binaries can be destroyed without af- 
fecting the shape of the mass ratio distribution. This subtle 
dependence of binary destruction on both the orbital and 
critical velocities for systems with unequal masses can be 
summarised by consid ering the hard-soft boundary a s the 
"fast-slow" boundary JHills fc Disslvll 19891 ; lHirlslll990r i. 



2 ANALYTICS OF BINARY DESTRUCTION 



Sever al dynamical classes of binary exist (jHeggial 19751 : iHilla 
ll975T ). 'Hard binaries' have a binding energy greater than the 
local Maxwellian energy ({m)a 2 , where (m) is the average 
stellar mass and a is the Maxwellian velocity), and are un- 
likely to be destroyed in an interaction (indeed, they are of- 
ten 'hardened'). 'Soft' binaries have a binding energy which 
is much lower than the local Maxwellia n energy, and can be 
disru pted by very small perturbations (jBinnev fc Tremainei 
Il987f ). Finally, 'intermediate' binaries have a binding energy 
comparable to the local Maxwellian energy. The binding en- 
ergy of a binary, -Ebind, can be written as 



Ebind = 



Gm p m s 



(1) 



where m p and m s axe the primary and secondary masses 
and a is the semi-major axis. If two binaries have the same 
separation, a, they should have slightly different binding en- 
ergies if their mass ratios are different. If we consider two 
binary systems with equal separations, but one is a 1 Mq- 
IMq system and the other is a IMq-O.IMq system, then 
we expect the latter system to have a lower binding energy 
(and hence may be more susceptible to destruction) than 
the former. 

However, analytic considerations of binary destruction 
are further co mplicated when the m ass ratio of the system is 
not unity (e.g. IFregeau et al.l l2006) . If we consider an inter- 
action between a binary and perturbing star of mass Mp ert , 
the critical velocity of the system (where the total energy of 
the binary and perturbing star is zero) can be written as: 



/ GnipT 



fia 



1/2 



(2) 



where fi is the reduced mass of the binary-perturbing star 
system: 



fJ 



(m p + m s )Mp CTt 



m p + m s + Mport 
Furthermore, if the orbital velocity, n or b is 

( G(m p + m s ) x ' '" 



Vorb = 






(3) 



(4) 



then for comparable masses (m p = m 3 — Mp crt ), v c ~ v OI b- 
However, if m p >> m s , but m p ~ Mp cr t, then« c ~ ^/2qv OY \y. 



3 METHOD 

3.1 Cluster structure and virial state 

We run four suites of simulations to examine the po- 
tential change in CMRD through dynamical interac- 
tions in a clustered environment. In the first three sets 
of simulations, the binary systems are placed in cool, 
clumpy clusters IjParker. Goodwin fc Allisonl 12011! ) because 
such initial conditions reflect the fact that many star 
forming regions are observed to be substructured (e.g. 
ICartwright fc Whitworthll2004l:ISanchez fc Alfaroll2009r i and 
subvirial (e.g. IPeretto et al. 2006; Fiir esz et al.l l2008). 

Furthermore, IParker fc Goodwinl i|2012r ) have recently 
shown that these cool, clumpy clusters can stochastically 
alter identical binary populations so that the separation 
distributions are statistically different after 1 Myr in differ- 
ent clusters. Therefore, these initial conditions represent the 
most extreme way of disrupting binaries and may alter the 
CMRD to an even greater extent than say, smooth, viri- 
alised clusters. We also run one suite of simulations with the 
binaries distributed in smooth, virialised Plummer spheres 
(Plu mmerlll911l ) as a control run. 

The substructured clusters are set up as fractals with a 
radius tf = 1 pc, according t o the prescription detailed in 
iGoodwin fc Whitworthl (120041 ). Fractals have the advantage 
that the amount of substructure is described by just one 
number (the fractal dimension, D), although it is unclear 
whether stars actually form in a fract al distribution (they 
may do: lElmegreen fc Emiegreenll200ll ). Highly fractal clus- 
ters (D — 1.6) process more binaries than smoother clusters 
e. g. D = 2.6) du e to the dense pockets of substructure 



Parker et al.l 120111 ). and in this work we adopt D — 2.0, 
which is a 'moderate' level of substructure. 

The stars in our fractal clusters have correlated veloc- 
ities, so that two nearby stars in the fractal have similar 
velocities, whereas dista nt stars may have very differen t ve- 
locities, as d e scribe d in IGoodwin fc Whitworthl |200J) and 
IParker et all (|201ll ). 

In the suite of simulations in which the stars are 
distributed in a smooth Plummer sphere, the positions 
and velocities are chosen accord ing to the prescription in 
lAarseth. Henon fc Wielenl l| 19741 ). 

Finally, we scale the velocities of the stars in the cluster 
to the desired virial ratio (Q v ir = T/|f2| where T and \Q\ are 
the total kinetic energy and total potential energy of the 
stars, respectively). For the fractal cluster simulations we 
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choose a subvirial (cool) virial ratio, Q v ir = 0.3, whereas the 
Plummer sphere clusters are in virial equilibrium (Q v ir = 
0.5). 



3.2 Initial binary populations 

We create the clusters so that every star is placed in a binary 
system i nitiall y. Primary masses are drawn from a 2-part 
iKroupal (120021 1 IMF of the form: 



N(M) oc 



M 
M 



-1.3 
-2.3 



m < M/Mq < mi . 
mi < M/Mq < m 2 . 



(•») 



where m = 0.1 M , mi = 0.5 M Q , and m 2 = 50 M . We 
choose secondary masses a ccording to the flat CMRD ob - 
served in the Galactic field JReggiani fc Mever|[2oTH . l2013h . 
For a comprehensive discussion of different binary compo- 
nent pairing algorithms, w e refer the interested reader to 
iKouwenhoven et al.l ((2009J). We do not allow brown dwarf 
primaries in the simulations. 

In two sets of si mulations we draw periods from the 
distribution f ound by Duquennov fc Mayor! (|1991f ) for field 
G-dwarfs fsee lRaghavan et aUl2010l . for an updated but sim- 
ilar fit): 



/ (logic, P) = exp 



-(logioP-logipP) 
2rr 2 



(6) 



where logioP = 4.8, ai ogl0 p = 2.3 and P is in days. We 
convert the periods to semi-major axes using the masses 
of the binary components. This corresponds to semi-major 
axes in the range 10~ 2 — 10 5 au. 

The eccentricities of the binary orbits are drawn from 
a thermal distribution (|Heggidl 19751 ) of the form 



Me) = 2e. 



(7) 



In the sample of lDuquennov fc Mavorl ql991f ). close binaries 
(with periods less than 10 days) are almost exclusively on 
tidally circularised orbits. We account for this by reselecting 
the eccentricity of a system if it exceeds the following period- 
dependent value: 



e tid = 77 [0.95 + tanh (0.6 logi P - 1.7)] 



(8) 



In the remaining two simulations, we assign all bina- 
ries a semi-major axis of either 30 au or 10 au, in order to 
test whether the shape of the CMRD could be altered by 
dynamical interactions when the semi-major axis becomes a 
constant in Equation [l] 

We combine the primary and secondary masses of the 
binaries with their semi-major axes and eccentricities to de- 
termine the relative velocity and radial components of the 
stars in each system. The binaries are then placed at the 
centre of mass and velocity for each system in either the 
fractal distribution or Plummer sphere (see Section [371}. 

The simulations are run for 10 Myr using the kira in- 
tegrator in the Starlab package (e.g. IPortegies Zwart et al.l 
ll999j, |200lj). We do not include stellar evolution in the sim- 
ulations. A summary of the four sets of simulations is given 
in Tabled] 



Table 1. A summary of the four sets of simulations. The values 
in the columns are: the number of stars in each cluster (Altars), 
the morphology of the cluster (either a Plummer sphere or frac- 
tal), the initial virial ratio of the cluster (Q v ir)i the initial ra- 
dius of the fractal, (rp), or the initial half-mass radius of the 
Plummer sphere, (ri/2), the fractal dimension, D (if applica- 
ble) and the adopted semi- major axis distribution, a-^- ln (cither 
iDuquennov fc Mavorl 111991!) , a delta function at 30 au or a delta 
function at 10 au). 
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4 RESULTS 

In this section we first compare the shape of the semi- 
major axis distribution and companion mass ratio distri- 
bution (CMRD) at and 10 Myr (i.e. before and after dy- 
namical evolution in the cluster). We then consider the in- 
teraction histories of binaries that are destroyed. 



4.1 Initial and final distributions 

We compare the initial (OMyr) semi- major axis (hereafter 
separation) distribution and CMRD to the distributions af- 
ter 10 Myr. Most binary destruction in the ty pe of clusters 
we simulate here occurs i n the first 1 Myr (e .g. IParker et al.l 
120091 : iMarks et ai]|201ll : IParker et all 1201 J ), but we anal- 
yse the simulations at 10 Myr to ensure that all destructive 
encounters are tracked. After 10 Myr the clusters are so dif- 
fuse that very little further dynamical processing will occur 
before the binaries become members of the Galactic field 
population. We determine whether a binary is energetically 
b ound or not us i ng the nea rest neighbour algorithm ou tlined 
in IParker etafl (|2009h and IKouwenhoven et ail (|2010T ). The 
number of binaries, and the binary fractions at and 10 Myr 
are given in Table [5] for each simulation. 

In our first simulation, we placed 750 binary systems 
with a field-like separation distribution in a cool, clumpy star 
cluster. Such clusters violently relax during their collapse to 
a centrally concentra ted, spherical cluster (see figs. 1 and 2 
in IParker et al.ll201ll ). and binaries are processed during the 
erasure of substructure and due to the dense core which 
forms. In Fig. [T] we show the evolution of (a) the separation 
distribution and (b) the CMRD. The initial distributions 
are shown by the red dashed lines and the final (10 Myr) 
distributions are shown by the solid black lines. 

Clearly, the shape of the separation distribution changes 
during the first 10 Myr of evolution. Binaries are destroyed 
(the binary fraction decreases from a primordial value of 
85 per cenu to 60 per cent after 10 Myr) , especially the 



1 Note that formally we placed every star in a binary initially. 
However, the initial cluster density causes the widest binaries to 
be unbound, even before dynamical evolution. These initially un- 
bound binaries are not included in the analysis. 
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(a) Separation distribution 



(b) Mass ratio distribution 



Figure 1. Evolution of the binary cu mulative separation d i stribu tion and cumulative companion mass ratio distribution (CMRD) for 
binaries with periods drawn from the iDuquennov &; Mavorl IJ1991T ) fit to the Galactic field data. The initial distributions at Myr are 
shown by the red dashed lines, and the distributions after 10 Myr of dynamical evolution are shown by the solid black lines. 

Table 2. A summary of the numbers of binaries destroyed in each simulation. The values in the columns are: the simulation number, 
the morphology of the cluste r (either a Plumme r sphere or fractal), the initial virial ratio of the cluster (Q v ir)i the adopted semi-major 
axis distribution, a^in (either IDuquennov fc Mayor (1991), a delta function at 30 au or a delta function at 10 au), the number of binaries 
at OMyr (-/Vbin,0Myr)i the binary fraction at OMyr (/bin, OMyr); the number of binaries at 10 Myr (iVbin.iOMyr) and the binary fraction 

at lOMyr (/bin, 10 Myr)- 
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654 
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512 
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wider systems (> 200 au) which are intermediate-soft in 
such dense environments according to the Heggie-Hills law 
|Heggielll975l ; lHillslll975l ) and are broken up. However, the 
CMRD retains the same shape, even after dynamical evo- 
lution. We find (almost) identical results in the simulations 
of Plummer-sphere clusters (comparison of the evolution of 
the binary fracti on and separation d istribution for Plummer 
sphere clusters ([Parker et al.l 2009) and for fractal clusters 
(jParker et al.ll201ll ) also shows that the amount of dynami- 
cal evolution in both morphologies is similar). 



In Fig. [2] we investigate whether placing binaries in clus- 
ters with identical (30 au) separations results in different 
initial and final CMRDs. Again, a certain number of bi- 
naries are broken up (the binary fraction decreases from 
100 per cent to 80 per cent after 10 Myr), and a smaller frac- 
tion have their separation hardened or softened (see panel 
(a) of Fig. [2]), but the shape of the CMRD is unchanged 
(Fig. [2(b)| . 



4.2 Binary disruption energetics 

These results indicate that destruction of binaries is inde- 
pendent of the mass ratio, q. However, according to Equa- 
tion [T] if two binaries have the same separation, they should 
have slightly different binding energies if their mass ratios 
are different. 

We demonstra te this in Fig. [3] wher e we d raw binaries 
randomly from the IDuquennov fc Mavorl (|199ll ) semi-major 
axis distribution, choosing the primary from the IMF and 
assigning th e secondary from a flat companion mass ratio 
distribution (Reggi ani fc Meverll201ll ). The solid line shows 
the cumulative distribution of binding energies for all bina- 
ries, and the dashed lines show the distributions for systems 
with q ^ 0.5 (the left-hand line) and q > 0.5 (the right hand 
line). 

The binding energy is weakly dependent on q, but of far 
more importance is the energy of a star which may break up 
the binary. A binary tends to be broken up in an encounter 
if the kinetic energy of the perturber exceeds the binding 
energy of the binary. We show the kinetic energy of a 1 Mq 
star travelling at (a) 1 kms -1 (the lefthand dotted red line) 
and (b) 20kms~ 1 (the right hand dotted blue line). Stars 
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Figure 2. Evolution of the binary cumulative separation distribution and cumulative companion mass ratio distribution (CMRD) for 
binaries with initial separations of 30 au. The initial distributions at Myr are shown by the red dashed lines, and the distributions after 
lOMyr of dynamical evolution are shown by the solid black lines. 
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determine the kinetic energy, J5 por t, of the interaction which 
destroys the binary: 



1 31 1 32 1 33 1 34 1 35 1 36 1 37 1 38 1 39 1 40 1 41 1 o 42 
Energy (J) 

Figure 3. Cumulative distribution of binary bindin g ener- 
gies f or systems with primary masses drawn from a iKroupal 
q2002l ') IMF, se condary masses drawn from a flat mass ra- 
tio distribution IIReggiani fc Meverl 120111) and periods from the 
iDuquennov fc Mavorl J1991T) distribution observed in the Galac- 
tic field. We show the distribution for all systems by the solid line, 
systems with q ^ 0.5 by the left-hand dashed line, and systems 
with q > 0.5 by the right-hand dashed line (as binding energy 
weakly depends on the mass ratio). We show the kinetic energy 
of a 1 Mq star travelling at (a) lkms' 1 with the lefthand dotted 
red line and (b) 20kms~ 1 with the right hand dotted blue line. 



in clusters typically have velocities ~lkms^\ but can suf- 
fer dynamical interactions which significantly increase their 
kinetic energy. Fig. [3] shows that the differences in binding 
energy due to low mass ratios are small compared to the 
range of kinetic energies a star can have in a clustered envi- 
ronment. 

With this in mind, we examine the distribution of en- 
counter energies that break up binaries. We track the inter- 
action history of every binary in our simulations, and then 



1 2 

^pert = 7T P crt ^ , P Grt " 



(9) 



Here, M port is the mass of the perturbing star, and v paTt is 
the magnitude of the relative velocity of the perturbing star 
with respect to the binary. If the perturbing star happens to 
also be a binary, we calculate v pclt using the centre of mass 
velocity of that binary, and use the mass of the binary for 
M port • 

In Fig. [4] we show the distribution of the kinetic energy 
of the interaction between a perturbing star and a binary 
which destroys the binary, divided by the binding energy 
of the binary. The open histrogram is the distribution for 
the simulation which has b inary separations drawn from 
IDuquennov fc Mavorl (|l99ll ). the hashed histogram is for 
the simulation which has all binaries with 30 au separations 
and the solid histogram is the simulation where all binaries 
have 10 au separations. The histograms are normalised to 
the number of binaries that are destroyed between and 
10 Myr in the first simulation (with the field-like separation 
distribution) to demonstrate the decrease in the number of 
systems that are destroyed with decreasing binary separa- 
tion. 

For the simula ti on w ith separations drawn from 
IDuquennov fc Mavorl (|l99ll ). the distribution of perturb- 
ing kinetic energy to binding energy peaks at a ratio 
Spert/Sbind ~ 10 — 100, suggesting that the input energy 
required to destroy a binary typically exceeds the binding 
energy by a factor similar to, or in excess of the most ex- 
treme mass ratios. However, the distributions for the simu- 
lations with separations drawn from delta functions of 30 au 
and 10 au (the hashed and solid histograms, respectively) 
indicate that for lower separations, the ratio of E pel t/ E^ind 
decreases towards unity. This could in principle lead to evo- 
lution of the CMRD, as the peak of the distribution of 
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Figure 4. Distribution of the kinetic energy of the interaction be- 
tween a perturbing star and a binary which destroys the binary, 
divided by the binding energy of the binary, for three of our simu- 
lations. The open histrogram is the distribu tion for the simulation 
which has binary separations drawn from iDuquennov fc Mavorl 
1119911 ) . the hashed histogram is for the simulation which has all 
binaries with 30 au separations and the solid histogram is the sim- 
ulation where all binaries have 10 au separations. The histograms 
are normalised to the number of binaries that are destroyed after 
10 Myr in the first simulation. 



-Epert/Sbind moves to values similar to the difference in bind- 
ing energy between systems with vastly different mass ratios. 
In practice, this is unlikely; so few binaries are destroyed in 
the cluster with binary separations at 10 au (typically only 
50 out of an initial total of 750) that a drastic bias towards 
destroying systems with low q would be required before this 
became evident in the distribution. 

A small fraction (usually less than 10 per cent, depend- 
ing on the initial separation distribution) of the destroyed 
binaries have an -Epert/Sbind ratio less than unity. These bi- 
naries are destroyed by multiple interactions, which have the 
progressive effect of lowering the binding energy and leav- 
ing the binary more susceptible to destruction after each 
encounter. 

Finally, we note that the shape of the CMRD could 
be altered by exchange interactions (for example, if up cr t < 
v c ). In this scenario, the least massive star in the binary is 
replaced by a more massive intruder, and the mass ratio will 
therefore increase. In our simulations, we find that typically 
< 10 exchange interactions occur per cluster, which is too 
small a rate to affect the shape of the CMRD. 



companion mass ratio distribution (CMRD) is unaffected by 
dynamical evolution in dense clusters. 

The (lack of) evolution of the CMRD is also unaffected 
by the assumed initial orbital parameters of binary systems. 
In the simulations presented in Fig. [2] we created binaries 
with separations drawn from a delta function at 30 au. Al- 
though we would expect the binding energy of binary sys- 
tems to be weakly dependent on the binary mass ratio, in 
practice the shape of the CMRD remains constant. The rea- 
sons for this are twofold. Firstly, the average energy of per- 
turbing stars (travelling on average at ~lkms _1 ) in the 
clustered environment is higher than the binding energy of 
half of the binaries in the cluster, with only a weak depen- 
dence on mass ratio. Secondly, the interactions that do de- 
stroy binary systems typically have a perturbing energy of 
order 10 - 100 times the binding energy; i.e. they are so en- 
ergetic that the dependence of the binding energy on mass 
ratio is irrelevant. Finally, a binary may suffer multiple inter- 
actions in a clustered environment, and can be de-stablised 
by previous interactions, before the encounter that finally 
destroys it. 

As we have shown in Fig. 2] the harder a binary, the 
lower the ratio of perturbing energy to binding energy is. 
Therefore, there may be a binary separation regime (< 10 au 
in the clusters we simulate here) in which the CMRD could 
change. However, at such low separations very few bina- 
ries are destroyed, and so the shape of the CMRD would 
be expected to remain roughly constant. A similar paucity 
of exchange interactions in these clusters also prevents the 
shape of the CMRD changing towards larger q values. 

Previous work on dynamical evolution of binaries in star 
clusters suggested that the shape of the CMRD is altered by 
dynamics. Kroupa ( 1995a) plotted the initial and final dis- 
tributions of the secondary mass (m s ) component of binaries 
with a G-type primary (0.8 < rrip/M© < 1.2) and showed 
that more binaries with m„ ~ 0.2 Mq are destroyed overall 
(see his fig. 4). However. iKroupal (|1995aT ) chose secondary 
masses from randomly sampling an IMF, the peak of which 
lies around (m) = 0.2 M©. Therefore, it is unsurprising that 
most binaries that are destroyed have a relatively low-mass 
companion and such a plot does not show the evolution of 
the mass ratio distribution. When one compares the shape 
of the CMRD, however, it is clear that dynamics does not 
alter an initial distribution. 

More recent work by iFregeau et al.l (J2004D examined 
binary-binary and binary-single star interactions in glob- 
ular clusters, and found that destruction is only very 
weakly dependent on the binary mass ratioqfj. In di- 
rect iV-body models of the old open cluster NGC 188, 



5 DISCUSSION 

Using iV-body simulations, we follow the dynamical evolu- 
tion of clusters containing 750 primordial binary stars for 
10 Myr. In our simulations of clusters, which are initially 
subvirial and substructured, binaries are destroyed imme- 
diately on local scales (i.e. within the substructure), be- 
fore the global collapse of the cluster at around ~0.8Myr 
fe.g. lAllison et alJl20irj|Parker et alfcoill ). The results pre- 
sented in Section|4]indicate that the shape of the binary star 



2 A wealth of literature also exists which presents numerical scat- 
tering experiments examining the effects of binary-binary and 
binary— single star interactions (e.g. iHut fc Bahca ll 1983; Mikkola 
1 19831 ; lHillsill990l ; iHeggie fc Hulll993h . However, these experi- 
ments tend to focus on systems with equal mass (m 



1u 
s 



1) binary components, though see ISigurdsson fc Phinnevl 
il993h ; IHeggie, Hut fc McMillan] dl996h for a description of ex- 
change interactions with unequal mass components. This, com- 
bined with a lack of information on subsequent encounters (which 
may further de-stablise soft binaries and could be important in a 
clustered environments) makes it difficult to compare the results 
to JV-body simulations of binary destruction in clusters. 
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iGeller. Hurley fc Mathieul (|2013T ) also found the initial and 
final CMRDs for main sequence binaries to be statistically 
indistinguishable. 

It should also be noted that the clustered environments 
modelled here have high densities (by design, so that some 
binaries are destroyed and we can examine the change in 
the C MRD) compared to t hose in the local Solar neighbour- 
hood (jBressert et alj|2010n , where the Orion Nebula Cluster 
i ONC ) is th e most dense region. Bas ed on the lBressert et all 
20101 ') data. lParker fc Meveri < |2012f> recently estimated that 
up to 50 per cent of star forming events could be dense 
enough to affect binary systems, so it is possible that many 
binary systems do not undergo significant dynamical pro- 
cessing. We have shown here that even in clusters with sim- 
ilar densities to the ONC (extrema in terms of local star 
formation), dynamical interactions are unlikely to affect the 
shape of the CMRD. 

Furthermore, we have not applied any cuts to our sim- 
ulations in terms of primary mass, separation or mass ratio 
ranges, which are often applied to observationa l samples to 
make them consistent (IReggiani fc Meveri 1201 If ). Typically, 
observatio ns of visual binarie s in clusters span the range 10s 
- 100s au (|King et alj|2012al lbl). and we have already seen 
that the shape of the CMRD is not affected by dynamical 
evolution in this separation range. If we were to impose a 
primary mass range cut, this would not affect our results as 
such a cut would bias our sampled q toward unity, thereby 
negating the argument of systems with lower binding ener- 
gies being more susceptible to destruction. 

There is observational evidence for a universal CMRD 
over a wide range of prim ary masses and q values 
(jMetchev fc Hillenbrandll2009l ). This appears to be true both 
in the field and in some associat ions and star forming regions 
JReggiani fc Meveri [20TH . 120131 '). However, due to its inde- 
pendence of dynamical evolution, the CMRD may help trace 
differences between star formation events and test models of 
binary formation. If differences do exist, the combined study 
of the CMRD in different environments can be also used to 
determine which types of clustered co nfiguration con tribute 
most to the field population (see also lGoodwinll2013h . 

Finally, we note that uncovering evidence for a change 
in the shape of the field CMRD towards random pair- 
ing for wide binaries would support the hypothesis of 
wide bin ary formation during the dissolution phase of star 
clusters JKouwenhoven et all I20T0I ; iMoeckel fc Bate! l20ld ; 
iMoeckel fc Clarke! l201lh , as such systems effectively form 
via capture rather than core or disc fragmentation. 



energy than systems with higher mass ratios, and hence be 
more susceptible to destruction. This is not a strong effect, 
however, and binaries are typically destroyed by a perturber 
with kinetic energy well in excess of the binding energy. 

(iii) Any differences in the observed CMRD between dif- 
ferent star forming regions would indicate different modes of 
star formation, and hence whether star formation is univer- 
sal or not. The CMRD is therefore a stronger diagnostic for 
searching for different modes of star formation than using 
e.g. the binary separation distribution, because the shape of 
the separation distribution does change through dynamical 
interactions. 
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